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The present paper is a review of the phenomena related to non-equilibrium electron relaxation 
in bulk and nano-scale metallic samples. The workable Two- Temperature Model (TTM) based on 
Boltzmann-Bloch-Peierls (BBP) kinetic equation has been applied to study the ultra- fast (femto- 
second) electronic relaxation in various metallic systems. The advent of new ultra-fast (femto- 
second) laser technology and pump-probe spectroscopy has produced wealth of new results for 
micro and nano-scale electronic technology. The aim of this paper is to clarify the TTM, conditions 
of its validity and non-validity, its modifications for nano-systems, to sum-up the progress, and to 
point out open problems in this field. We also give a phenomenological integro-differential equation 
for the kinetics of non-degenerate electrons that goes beyond the TTM. 

PACS numbers: 63.20.Kr,72.10.Di,72.15.Lh,72.20.Dp 

"The first processes, therefore, in the effectual studies of the sciences, must be ones of simplification and 
reduction of the results of previous investigations to a form in which the mind can grasp them." — J.C. MAXWELL 

I. INTRODUCTION AND BRIEF HISTORY 

The Two- Temperature Model (TTM) describes a non-equilibrium state between electrons and the lattice. When an 
electric field is applied to a metallic sample or the sample is photoexcited by a laser pulse, a state of non-equilibrium 
arise between electrons and the lattice. Since the time required to establish equilibrium in the electron gas (degenerate 
Fermi distribution) is much less than the time required to establish equilibrium between the electrons and the lattice 
phonons, the metal can be considered as composed of interacting subsystems, one electrons and other phonons. The hot 
degenerate (thermalized) electron gas relaxes to the bath phonons with relatively slow electron-phonon interactions. 
This problem was solved by V. L. Ginzburg and V. P. Shabanskii (1955) 1 in high temperature limit (T 3> To, where 
Tjj is the Debye temperature). The complete solution for arbitrary temperatures was given by Kaganov, Lifshitz 
and Tanatarov (1956) 2 . They considered the weakly coupled electron and phonon subsystems and ignored the band 
structure of the metal, this work was further extended to the case of a metal exposed to ultrashort laser pulses^, 
and today known as TTM. It is a pioneering work, that provides the basic model for ultra-fast pump-probe laser 
studies in metals^' 5 -' 6 - ' 7 ' 8 ' 12 ' 13 . The theory of thermal relaxation of electrons in metals was further extended by P. B. 
Allen 14 . It has been proposed that TTM can be derived from D. N. Zubarev's non-equilibrium statistical operator 
(NSO) metho d 15 ' 16 by systematic approximations, with this, one can obtain even more advanced results than TTM 
including the non-Markovian effects^. NSO method is based on a generalization of Gibbs' method of equilibrium 
states to non-equilibrium states. Gibbs equilibrium statistical mechanics is based upon the deep connection between 
additive integrals of motion (most common is energy) and the invariance of phase space distribution function due to 
Liouvelle theorem. This idea has been generalized to the non-equilibrium states by D. N. Zubarev using the concept 
of local quasi- integrals of motion and N. N. Bogolyubov's idea of a hierarchy of relaxation times i.e., evolution of the 
system from far-from-equilibrium stage to kinetic stage and then to the hydrodynamic stage as time proceeds. In 
our context of non-equilibrium electron-phonon system these time scales are in pico-second domain. It is a hope that 
some one among the readers will accomplish this very ambitious task of obtaining TTM from NSO method. 

For the case of bulk semiconductors, theory of electron relaxation was given by Sh. M. Kogan (1963}i£, with the 
similar assumptions as in the case of metals. It was extended and refined by S. Das Sarma (1990) 20 and M.W.C. 
Dharma-Wardana (1991) 21 using non-equilibrium Green's functions with dynamical screening and hot phonon effect. 
Large amount of work has been done and going on in nano-scale materials, the advent of nano-technology and ultra- 
fast lasers has completely revolutionized the field. However, this review is devoted to metallic nano-systems only. In 
metallic nano-systems the work of P. M. Tomchuk 2 ^ is quite important (see the review^ 3 , for more detail on various 
phenomena caused by hot electrons in island metal films). 

The aim of this paper is to clarify the TTM, conditions of its validity and non- validity, and modified TTM as applied 
to nanoscale metallic systems, starting from first principles, without phcnomcnologically introducing power transferred 
to surface phenons as done in previous worksS 2 -. The second aim of this paper is to point out new open problems about 
non-equilibrium electron relaxation in metallic nano-systems(see last paragraph of section III B), and also in bulk 
systems (beyond TTM, section IV). We also give a simple stochastic model to describe relaxation phenomena in non- 
equilibrium, non-degenerate (non-boltzmann) electronic subsystem. The adjective "non-Boltzmann" signifies that the 
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TABLE I: Notation used for important physical quantities 
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T^vnicftl valup ftnH uriitt'^ 


Temperature of the hot thermalized electron distribution 




~ j Ann is 
~ 4UU rv 


Temperature of the phonon bath 


T 


~ 300 K 


Electronic heat capacity 


/ — 1 rri 


7 = bo Jm K 


Phononic heat capacity 


Cp 


~ 100 C e 


Electron-phonon interaction constant 


u b 


~ 10~ 19 Joule 


Electron-phonon interaction coefficient 


G 


~ 5 x 10 16 Jm- :i Sec- 1 K- 1 


Decay time of non-thermal electron distribution 


' e— e 


~ 300 femto-seconds (fs) 


Rise time of thermalized electron distribution 


Tth 


~ 500 fs 


Decay time of thermalized electron distribution 


7~e—p 


~ 1 — 2 peco-seconds (ps) 



gas is in non-equilibrium condition even in the classical approximation. We emphasis both aspects (namely, modified 
TTM for nano-systems using first principles and beyond TTM i.e., a general stochastic model) of the relaxation 
problem. 

The paper is planned as follows. In the next section the approximation of metal as a TTM system is introduced and 
clarified. Section III is devoted to electronic relaxation in metallic nano-systems where TTM is used with appropriate 
boundary constraints. In section IV, by pointing out the limitations of TTM, we introduce a generalized stochastic 
model for electron-phonon interaction. In section V, we will end the manuscript with brief conclusion. 

Before we go into the details of TTM, we first carefully define our notation. There is a lot of notational ambiguity 
in this research field, notation varies from author to author. However, for the sake of a coherent presentation we fix 
our notation as tabulated in table I, also, units and typical values are given for a physical feeling about the quantities 
involved. The terms non-degenerate, non-Fermi-Dirac, or non-thermalized electron distribution has been used in the 
literature, all have the same meaning, and point to the same physical condition of the electrons, but here, we adopt 
the term "non-thermalized electron distribution" for the sake of a unified presentation. 

II. WHAT IS TWO-TEMPERATURE MODEL(TTM)? 

Consider a metallic sample photoexcited by a femto-second laser pulse, because of the large difference between the 
electronic (C e ) and lattice(phononic C p ) heat capacities (with C p ^> C e at room temperature), the femto-second laser 
pulse creates non-equilibrium electron distribution, leaving the lattice temperature essentially unchanged T ~ ZOOK . 
Then, over a time scale of hundreds of femto-seconds, the non-equilibrium electrons redistribute their energies among 
themselves through electron-electron coulombic interaction, and return to a local equilibrium (among themselves) at 
a somewhat elevated temperature T e > T. It is called the thermalized electron distribution. This excited thermalized 
electron gas then cools(relaxes) via the electron-phonon interactions, giving up the excess energy to the phonon bath. 
Thus, the widely separated time-scales(the intra-electron and the intra-phonon relaxation times <SC the inter-electron- 
phonon time scale) justifies defining the two temperatures T e and T. This motivates the Two- Temperature model 2 -! 3 -. 
The TTM describe this relaxation process, and has been used extensively by the workers in the field of ultra-fast laser 
spectroscopy in nano-scale materials 32 . Briefly, the TTM assumes 

(a) The electron-electron(coulombic) and the phonon-phonon(anharmonic) processes are much faster than the 
electron-phonon processes, so as to define T e and T(=/= T e in general) and maintain their local equilibrium distributions 
giving 

1 1 

(for electrons) Nu = ^ , e = - — — 

(for phonons) Nt — -37 ,8 — , (1) 

with a fermionic electron distribution at temperature T e and a bosonic phonon distribution at temperature T (T < T e ). 

(b) Homogeneous excitation and no spatial diffusion. 

(c) Delta-pulse laser excitation. 

So, only the collsion term is important in the Boltzmann transport equation, which has two contributions; phonon 
generation or phonon absorption when electron scatters from one state to another. For a general discussion of the 
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applicability of the Boltzmann kinetic equation in quantum cases see the review article by A. L. Kuzemsky 1 - 8 . The 
transition rate (probability/second) of an electron scattering from a state k to state k' with a phonon generation of 
energy hvf is W(k'\k)S(e k ' — e k — and for phonon absorption, transition rate is W(k\k')S(e k — e k > + hu>). 



W(k\k') = 



P VS 2 



■Uf 



f = k-k'. 



(2) 



With Ub as the electron-phonon interaction constant (it appreas in the expression of time of flight of electrons and 
is a measure of interaction energy). It should be clearly distingueshed from electron-phonon interaction coefficient 
G (see table I), p, V and S is the metal density, unit cell volume and sound speed respectively. So, the net phonon 
generation rate per unit volume is 



N 



f 



W(k\k'){(N f + 1)N„(1 - N k ) - N f N k (l - N k ,)}S{e k , -e k - hu)(2 / (2vf)dT k , . 



(3) 



This equation is known as Boltzmann-Bloch-Peierls (BBP) kinetic equation^ 4 -. The delta function ensures the energy 
conservation. The factor of 2 near the volume element is for electron spin degeneracy. Here dr denotes k-space volume 
element, but later on t is used to represent relaxation times, but confusion should not arise. One can show that the 
rate of energy transfer per unit volume by the electrons to bulk phonons is 2 , 



bulk 



dr -f 

NfhLUfV—4^ 
1 1 (2tt) 3 



m 2 [/ 2 w 4 £;F 



2(2ir) 3 h i pS 4 



[T e — T] . 



This can be cast in the following form 3 



d{C e T e ) 

dt 
d(C p T) 

dt 



= —G(T e - T), 
= G(T e — T). 



(4) 

(5) 
(6) 



The above coupled differential equations are the defining equations of the Two- Temperature Model (TTM) of hot 
electron cooling. G is the electron-phonon interaction coefficient. If the system is not delta-pulse excited and heating 
of electron gas is going on in parallel with relaxation processes, then one can add a term Q to the right hand side of 
equation (5) corresponding to the specific power absorbed by the sample from the laser field. 



III. NANO-SYSTEMS: TTM AND BOUNDARY CONSTRAINTS 



Nano-sized metallic systems are quite important for micro- and nano- electronic technology. In a typical electronic 
circuit the passage of current creates a situation of non-equilibrium between electrons and the phonons. In bulk 
metals, the electron energy relaxation is due to the cherenkov generation of acoustic waves. But this mechanism is 
not present in case of nano-sized metallic systems- 24 . TTM can not be directly applied to such small systems, because 
of the non-resonant nature of interaction between bulk phonons and the electrons. In nano-systems the main channel 
of electron energy transfer is the electron surface-phonon interaction(cffective electron mean-free path becomes of the 
order of particle size). In the following we assume that TTM holds good with extra geometric constraints due to 
nano-size. 



A. Nano-scale Metallic Films 



Here we extend the TTM to nano-scale metallic films by considering surface phonon generation. In nano-scale 
metallic films the electron-surface-phonon interaction is important because film thickness is of the order of electron 
mean free path. The state of electrons and phonons is described by equilibrium Fermi and Bose functions with 
different temperatures. The new feature one has to consider is the geometric constraint of surface phonons. We 
obtain expressions for the energy transfer rate from thermalized hot electrons to surface-phonons, which is order of 
magnitude less than that for the bulk. The whole process occurs at pico-second time scales. We consider the case of a 
homogeneously photoexcited nano-scale metal film, and strong damping of surface phonons, due to the strong coupling 
with the substrate on which film was developed. Thus, there is no surface standing modes at the film surface^. 

Consider a hot Fermi electron distribution at temperature T e and a surface phonon(2-D) equilibrium distribution at 
temperature T (T < T e ). Both thermalized distibutions are weakly interacting subsystems. The energy flows from hot 
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degenerate electron distribution to phonon bath. In the following, we calculate, in line with Kaganov etal^the energy 
transfered per second per unit volume from the hot thermalizcd electron distribution to the relatively cold phonon 
gas. The equilibrium distributions for electrons and phonons are given by equation (1). Energy and momentum 
conservation conditions gives, 



£fc' - £k = , k' k = f ,k' x -k x = f x , 
ky ~ ky ~ fy i = ~k-z j w = sj , 
= *s£, £(2[fc' x / x + fc' y / y ]-/ 2 )=^/, 



(7) 



On simplifying, 



2m 



[2fc'sin6»cos(> - 0') - /] = hs. 



(8) 



Where is the angle between fc^-axis and plane of incidence. <f>' is the angle between scattered phonon direction and 
fc^-axis, and 9 between incident electron direction and k z direction as shown in Figure 1. The probability W per unit 
time that the electron in a state with wave vector k' will scatter to a state with wave vector k by emitting a phonon 
of wave vector / is; 



W(k' - /; k') = aS(s k , -e k - Jiu), a = (nil* / pV S 2 S ) . 



(9) 



With U s as the electron surface-phonon interaction constant. Here p, V and S s is the metal density, unit cell volume 
and surface sound speed respectively. The change per unit time per unit volume in the number of surface-phonons 
with wave vector / and energy fiuj is (Boltzmann-Bloch-Peierls kinetic equation); 



N f = J au f {{N f + l)N k ,(l - N k ) - N f N k {l - N k ,)}6{e k , - e k - Tilo)(2/ (2Tif)dT k , 
Using the energy and momentum conservation equations, the delta function can be written as 



2m 

S{e k , - e k - M = - 5 - < y[(2fe / sin 9 cos(0 - 0') - /) 
ft / 



2ms 
~7p' 



and for a metal, we have (/ ~ 10 9 m _1 ) > ~ 10 7 m _1 ). So the equation (10) is 

N f = [au f {(Nf + l)N k ,(l - N k ) ~ N f N k (l - N k ,)}5(2k' sin6cos((f> - <jf) - /) 



4:772 fc' 

x — - sin 9d9d(f>dk' . 



27r 2 r/ 



(10) 



(ii) 



(12) 
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The last integral in the above equation is 

r-2n 



[(N f + l)JV fc ,(l - N k ) - N f N k (l - N k ,)]k' 2 dk' 



2tt 



S[2k' sin6»cos(^- <j>') - /]#. 







S[2k' sin0cos(0 - 0') - f]d<f> = — — - 



1 



|sin(^i - <£')| Jo 
1 



2tt 



fc'sin6>VTl - / 2 /(4/c' 2 sin 2 6»))' 
inserting this in (13) we get 

mas 



N 



f 



(2-Khy 



2tv 



6(fc - 0')rf0 



[(N f + l)N k ,(l - N k ) - N f N k (l - N k ,)]k'dk' 



(13) 



(14) 



(15) 
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The above mentioned process will always happen as from energy and momentum conservation, sin# ~ f/2k', which 
holds good in a metal as / < kl . By inserting for N e and Nf in equation (10) we get 



N f = 



mas 



(2-Khy 



1 



e B {E' k -huj f -e o ) k i dk t 



(16) 



Here, we will make an approximation to solve the integral in the above equation. The first approximation is that the 
phonon energy ?kj/(~ meV) <C fc_Br e (~ eV), the electron energy. So, /3 e Ku>f ~ 0. With this, the integral in Eq.(16) 
is 



1 



e /(fcm) + 1 



f{k m ) — Pe 



2m 



1m 



As 

mui 



■C |/9eo|i the quantity in the square brackets is order of unity. Finally, the integral in Eq.(16) is 



~ ^(l/( e ^ ft <"/ - 1)). With all this, Eq.(16) takes the form 



N 



f 



m asujf 



{2irh) 2 h 



( e Phu f _ i^ e e hw; - l)y ' 



(17) 



Here a is defined in equation (9). The energy transfered by the electrons to the surface-phonons per unit volume per 
unit time is 



U, 



surface 



(2^) 2 Jo 



Id, 



N f hu; f 2irfdf, 



(18) 



where fo s and a are the Debye wave vector for the surface phonons and lattice constant respectively. From Eq.(17) 
and Eq.(18) with relations u>ds = S s fos , Tt-^Ds = ksTus and setting x — fuOf /kBT e , we get 



U, 



surface 



nU 2 m 2 



(2Tr) 3 h'paSI 

4 r T Ds /T s l 



-dx 



T 



4 r T Ds /T x 3 



-dx 



(19) 



Here, Tjj s is the surface Debye temperature. Equation (19) can be simplified in two special cases, first, for low electron 
and phonon temperatures as compared to Debye temperature, i.e., T,T e <C Tr> s , Eq.(19) reduce to 
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{2nfh'paSl 



kBT Ds 
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T„ - T 4 
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-dx. 



1 



(20) 



An important point to be noted in the above equation is that the electron to phonon energy transfer 
rate depends upon 4 th power of electron and phonon temperatures as compared to the corresponding 
case in the bulk(there it is 5 th power of electron and phonon temperatures 2 ). In second special case, when 
T e ,T^>T Ds , we get 



surface 



TrU 2 m 2 



3(27t) 3 ?T paSl 



( k B T Ds 
{ U 



4 r 



T P -T 



Tds 



(21) 



The above equation (Eq.(21)) is the basics of what is called the TTM. The surface Debye temperature Tjj s — -r-fos, 
for two acoustic modes per atom is given by 



(27T) 2 Jo 



2nfdf — 2A surface, 



which gives fos — \/87m 2 / 3 , n is the number density per unit volume. Now, for the bulk case^ 



Ubulk = 



2(2Tr) 3 h 3 pS b l 



[T e -T]. 



(22) 



(23) 
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FIG. 1: Scattering of an electron from a surface. 



From Eq.(21) and Eq.(23) we have 




ace 




U s ] 2 nSi 



(24) 



U b \ auj% b 



For a gold metal film, assuming U s — Ub, with a = 4.1 x 10 _10 m , p = 19.3 x I0 3 kg/m 3 , n = 5.9 x 10 28 m~ 3 , To — 
185 K , lud = 2.42 x 10 13 rads / sec, the above ratio is ~ 0.088 or about 9 percent. 

In brief, we have treated here the problem of energy relaxation of photo excited degenerate electrons in a nano-scale 
metal film. The new feature one has to add is the geometric constraint of surface phonons which is responsible for 
reduced energy transfer rate in nano-scale metal films as compared to bulk metals. In this simple calculation, the 
effect of electron-phonon screening and quantum confinement effects due to nano size, are not taken into account. 



As explained in the previous section, TTM can not be directly applied to such small systems, because of the 
non-resonant nature of interaction between bulk phonons and the electrons. The main channel of electron energy 
transfer is the electron surface-phonon interaction(effective electron mean-free path becomes of the order of particle 
size)21 P. M. Tomchuk and E. D. Belotskii give a quantum-kinetic treatment of hot electron energy relaxation in 
small metallic particles^. They phenomenologically introduced the power transferred from hot electrons to surface 
phonons (acoustic and capillary modes on the surface) , and relate that to the microscopic parameters of the electron 
gas. In the following, we will give a calculation starting with BBP kinetic equation with geometric constraints, and 
obtain the expressions for electron surface-phonon coupling coefficient at both temperature limits. 

The reduced dimensionality of nano-systems give rise to two important physical effects. One is quantum size effect, 
which give the quasi-continuous energy spectrum, and another is the geometric size effect (surface area/volume) which 
enhances the electron-surface interaction and heat diffusion to the bath. So, the following points are important: 

(1) in quantum size regime when the particle size is typically less than 5nm, the band structure splits into discrete 
levels and the equilibrium partition function of electrons will not be the same as that for the bulk. The function 
depends upon evenness or oddness of the number of electrons in the particle^; 

(2) Since electronic mean free path(several hundred angstroms in metals) is more than the particle size, even at 
high temperatures, the scattering events from the surface of the particle will take place. If the time between two 
scattering events is less than the electronic internal thermalization time, one has to use non-equilibrium distribution 
functions to describe the problem of energy transfer from electrons to phonons^; 

(3) In quasi-continuum regime(particle size more than 10nm)the main channel of electron energy loss is through 
electron surface interaction. One can use the Two- Temperature Model for electrons and the surface phonons. The 
only extra conditions to be imposed are of geometric nature. The present calculation is done in quasi-continuum 
regime, considering point (3). The points (1) and (2) are not included in the present calculation. It is assumed that 
TTM holds good, but replacing bulk phonons by surface phonons. This will not be applicable for small time scales 
when the electron distribution function is not Fermi-Dirac. It is to be noted that, as the hot electrons lose their energy 
to the lattice, and after some time, the lattice will become very hot, this heating will reduce drastically the electron 
mean free path and cause the failure of the applicability of the model. But for the case of metal particles it takes 
about 2 pico-seconds to transfer the energy to the lattice bath, so the present model is applicable within this time 



B. Hot electron relaxation in a metal nanoparticle: electron surface-phonon interaction 
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scale. The dispersion relation used for the surface phonons is linear under Debye approximation and surface sound 
speed is determined in terms of elastic continuum theory assuming stress-free boundaries^. 

We again consider the case of a homogeneously(no spatial diffusion) photo-excited metal nanoparticle. The equi- 
librium distribution functions of electrons and phonons are defined by equation (1). 

We proceed on similar lines as in the previous section, and use the conservation of energy and momentum for one 
phonon scattering as shown in Figure 1. We use the BBP kinetic equation to calculate phonon generation rate. After 
imposing all boundary conditions we have 



N f 



aco f {{N f + 1)N«{1 - N k ) - N f N k (l - N k ,)}6(e k > - e k - huj){2V/{2^f)dT kl . 



(25) 



N f = 



m 2 asVojf 



{2irh) 2 h 



(26) 



Again a is defined in equation (9). Similarly, we calculate the energy transfered by electrons per unit volume per unit 
time to surface phonons. We use the elastic continuum theory assuming stress-free boundaries, for density of states 
on the particle surface^, Clearly 



U, 



surface 



N f hLu f (D s (Lu f )/((4/3)nR 3 ))du; f 



D s {w f ) 



cj f R 2 
2S 2 ' 



(27) 



Where u>ds is the Debye frequency for the surface phonons and D s {uS) is the surface phonon mode density^. Finally 
we obtain 



U, 



surface — I ^ f 





Vl = 



3am 2 V 



dujf 



_32n 3 h 2 S s R. 

The surface Debye frequency and the surface Debye temperature from 

D s {ujf)dujf = 2{AnR 2 a/V), 



(28) 



(29) 



are, lod s = y327raS| JV and Td s = (TiS s /kB)\/'i2'Ka/V respectively. With this we obtain 
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4 rT Ds /T x 3 



1 
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(30) 



As before, the above equation (30) can be simplified in two special cases, first, for low electron and phonon temperatures 
as compared to Debye temperature, i.e., T, T e <C Tjj s , we have 
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(31) 



Again, the point to be noted in the above equation is that the electron to phonon energy transfer rate depends upon 
4 th power of electron and phonon temperatures as compared to the corresponding case in the bulk(there it is 5 th 
power of electron and phonon temperatures^) . In second special case, when T e , T> Tds , we get 
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surface 



( k B T Ds \^ 


T e — T~ 


\ n J 





(32) 



The important electron surface-phonon coupling coefficient (see equation (5)) in femtosecond pump-probe experiments 
for nanoparticles is 



G 



3(V32^)m 2 U 2 k 
nh 3 pV 
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R' 



(33) 
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For a gold nanoparticle of radius R — lOnm, with U s — 10 19 joule , a = 4.1 x 10 10 m , p — 19.3 x I0 3 kg/m 3 , n = 
5.9 x 10 28 m- 3 , T D = 185 K , u) D = 2.42 x 10 13 rads/sec, 

G ~ 7.1 x 10 13 jowZe m^sec" 1 ^ 1 (34) 

Which agrees with experiments^ and P. M. Tomchuk's calculation. The electron-phonon coupling coefficient for the 
case of bulk is ~ 5 x 10 16 joule m~ 3 sec^ 1 K^ 1 . So, G (surface) is less by a factor of 10 3 from that of bulk, which 
indicates suppression of electron energy transfer to phonons in case of nano-particles. The present results show that 
the electron surface-phonon coupling constant will increase with the reduction of the particle size. So the hot electron 
thcrmalization time will reduce(fast relaxation) with decreasing size of the nanoparticl o 30 i 31 . The calculation does not 
include the effect of electron surface-phonon screening, but the fact that, due to electron wave function spill out and 
d- electron localization 3 ^ in nanoparticles, the screening will be comparatively less as compared with the bulk. The 
question regarding the weight of the two factors, namely, surface to volume ratio (geometric factor), and reduction of 
electron phonon screening, in the thermalization of hot electron distribution is still open. 



IV. BEYOND TWO-TEMPERATURE MODEL 



As discussed in section II, TTM assumes "instantaneous" thermalization of non-thermal electron distribution to 
a hot Fermi distribution, and then this hot Fermi distribution cools through electron-phonon interactions with the 
phonon subsystem. The assumption of "instantaneous thermalization" or "adiabatic assumption" is, in general, not 
in agreement with experimental observations^^, when one probe the system with very short and multi-wavelength 
laser pulses (~ 100/s). Thus the TTM is clearly inadequate to account for non-thermal electron distributions. First 
phcnomenological model to account for non-thermal electron distribution is given by Sun etal 8 , which is quite successful 
in explaining qualitatively the main features of the relaxation mechanism, and it is in agreement with experimental 
observations (typical pump-probe experiments) and numerical solutions of the Boltzmann equation. The model is 
applicable in the perturbative regime, where the system's response is linear and the measured changes in the reflectivity 
of the sample (thin Gold films of thickness ~ 100A) can be related to the changes in electron distribution. The main 
idea is that the transient reflection or transmission of the sample is proportional to the real and imaginary parts of 
the dielectric constant, which is further related to the changes in the electron distribution-through joint density of 
states-due to a model developed by Rosei etal*. Thus, in the linear perturbative regime, by knowing about transient 
reflectivity or transmissivity of the sample one can study the changes in the electron distribution. The rate equation 
model developed by them consists of dividing the whole sample into three interacting subsystems (see figure 2) (1) 
a very-low-density non-thcrmalized electron distribution, (2) high-density thcrmalized (Fermi) distribution, and (3) 
the phonon subsystem. Neglecting the particle exchange, the three subsystems interact according to: 

dN 

— = -X 1 N-X 2 N, 
at 

BT 

C e -^=-G(T e -T) + X 1 N, 

C p -^ = G(T e —T) + \ 2 N. (35) 

Here, N stands for the energy density stored in non-thermal part (subsystem 1), C e and C p are the electronic and 
lattice heat capacities, T e and T are corresponding temperatures, Ai is the energy loss rate from non-thermal part to 
the thermalized part of the electronic part and A2 = <j- is the loss rate to the phonon bath. Thus the energy flows 

N oc e~ l l r e-e from non-thermal part to the thermal part (subsystem 2) and also to the phonon bath (subsystem 3). 
In this process the temperature of the thermal part (subsystem 2) increases with time constant r t h and decays (due 
the electron-phonon interaction) with time constant r e _ p i.e., AT e oc (1— e _i / T *' l )e _ */ Te - t ', clearly, we have r t ^ > r e _ e 
(two loss channels from subsystem 1 and only one loss channel from subsystem 2). See tabic I for typical values of 
these time scales. 

In the perturbative regime (AT e <C Tq — equilibrium lattice temperature ~ T), the contribution of the thermalized 
electron distribution to the transient reflectivity of the film is proportional to the electron temperature change, and 
contribution of the non-thermal part is proportional to the instantaneous energy density N(t) of the non-thermal 
part. These two contributions to transient reflectivity depends upon the probe wavelength. This simple model yields 
remarkable agreement with the experimental results. However, a consistent micro-scopic model that take into account 
the general relationship between the optical response of the system and changes in electron distribution (not restricted 
in the perturbative regime) is not yet developed. The first principle model has to face the problem of non-equilibrium 
electron distributions, a part of the main problem of the non-equilibrium statistical mechanics itself. A possible way 
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FIG. 2: Excited metal consists of interacting subsystems. 



to the solution is to use D. N. Zubarev's non-equilibrium statistical operator method as mentioned in the introduction. 
But micro-scopic models exists for experiments with four- wave mixing of pulses where 'spatial parametric effect — is 
important and one have to consider non-Markovian effects^. 

To consider the non-equilibrium electron kinetics (non-thermal part), a stochastic kinetic model (SKM) is given 
below. By noting the fact that conduction electrons in metals can be treated free , an analytical treatment is 
given which is based on a generalization of the stochastic model known for a driven dissipative granular gas^. The 
generalized model is not applicable in an experiment where spatial parametric effect thus non-Markovian effects are 
important. The driven dissipative granular gas model is an interesting model where the particle-particle and the 
particle-bath collisions are parametrized in detail. More specifically , the total rate of collisions suffered by a given 
('tagged') particle is partitioned into the particle-bath collision rate (fraction /) and the particle-particle collision rate 
(fraction I — /). Further, a fraction a of the total energy of the colliding particles is partitioned randomly between 
the colliding particles, while the remaining fraction (I — a) is dissipated through the frictional contact during the 
collisions. The system is kept in the dynamic (non-Boltzmannian) non-equilibrium condition by a constant drive. 
In our generalization to the electronic system, the bath has the obvious identification with phonons, and the drive 
is to be identified with the photo-excitation. Also, the possibly dissipative electron-electron interaction has to be 
interpreted in terms of the coulomb interaction as screened by the dissipative polarization of the lattice. While, our 
generalization of the stochastic granular gas model to the electronic system covers time-dependent process relevant 
to the transient femtosecond photoexcitation, we have actually treated the steady state electron distribution under 
the cw(continuous wave) drive. 

This stochastic generalized model describes the scattering events for time less than r e _ e which could be quite large 
for the case of strong electron-phonon couplings. So, this model is quite general one. Also, for the case of continuous 
photoexcitation the total electronic system can be thought of divided into two subsystems, one non-degenerate and 
the other degenerate. This model captures the kinetic picture of the phenomenon in the non-degenerate subsystem. 



The stochastic kinetic model (SKM) for non-thermal electrons 

Let n e (E)dE be the number of electrons lying in the energy range ±dE/2 centred about E. The electron-electron 
collisions, assumed inelastic in general, are described by the process; Ei + E[ — ► Ef+E'j = a(Ei +E' i ) with a < I, in 
which the tagged electron of energy Ei collides with another electron of energy E[ lying in the energy shell E[ ± ^AE' i: 
and is scattered to the final state Ef. The scattering rate for this inelastic process is taken to be (1 — f )Tn(E' i )dE' i . 
Similarly, the electron-phonon scattering rate is given by fTripi^E'^dEl, with ripi^E'^dE^ as the number of thermal 
phonons in the phonon-energy shell E[ ± \AE[. Here, the fraction < / < 1 determines the relative strengths of the 
binary electron-electron and the electron-phonon collisions. Also, let the electrons are photoexcited at energy E ex 
above the Fermi energy at a rate g ex S(E — E ex ), and then recombine (deplete) from the non-degenerate distribution 
through recombination. This depletion rate can be modelled by a term ~gdS(E)n e (0). Here, Fermi energy is set 
equal to zero for simplicity. The phonons are assumed to remain in thermal equilibrium at temperatures T. For 
non-degenerate electronic subsystem energy to be the only label for the single particle states. The photo-excitation 
is taken to be homogenous over the sample, which is reasonable for a nanoscale metallic sample. For the above 
dissipative model driven far from equilibrium, the kinetics for the non-equilibrium electron number density n e (E) is 
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given by the rate equation 

= -n e {E) j dE'[n e (E')(l - /) + n ph (E')f}T 
+ J dzp(z) J dE' J dE"5{E - za{E' + E"))n e {E')n e {E"){\ - f)T 
+ J dz P (z) J dE' J dE"S(E - z(E' + E"))n e {E')n ph {E")fT 

+g ex (t)S(E -E ex )- g d S(E)n e (0). (36) 

In the above, we have assumed the total energy (E' + E") for a binary collision to be partitioned such that a 
fraction z, with probability density p(z), goes to the tagged electron of initial energy E' , and 1 — z to the colliding 
particle (electron or phonon of initial energy E"). The inclusion of a in the electron-electron collision takes care of 
the possibility of inelastic electron-electron collisions. Note that we have suppressed the time argument (t) in the 
non-equilibrium electron- number density n e (E). Taking the energy Laplace transform 

r-OG 

fie(«) = / e- sE n e (E)dE, (37) 
Jo 

of Eq.(36), we obtain, 

^n(s) = -Tn e (s)[(l-f)N e + fN ph ] + (l-f)r^ p(z)dz n 2 e (azs) 
+./T / dz p(z) n e (zs)n ph (zs) 



J o 

+g ex (t)e- sE °° -g d n e (0). (38) 

In the following, we will consider for simplicity the steady-state condition under constant (cw) photoexcitation, 
g e x(t) = g ex - A pulsed excitation can, of course, be considered in general. Accordingly, we set ^h e (s) = above, and 
all quantities on the R.H.S. of Eq.(38) are then independent of time. In order to calculate the steady-state electron 
distribution for the system in terms of the bath (phonon) temperature and other rate parameters, we expand h e (s) 
in powers of the Laplace-transform parameter s as 

n e (s) = N e - s{E e ) + s 2 (E 2 )/2..., (39) 

and equate the co-efficients of like powers of s. Thus, from the zeroth power of s, we obtain at once 

ne(0) = (g ex /g d ). (40) 

Similarly, from the first power of s, we get, 

(f/2)(e ph ) g e xE ex /T 

{ee) " Pe- Ph (i - «)(i - 75 + 772 + N$ hPe - ph [ Pe - ph (i -a)(i- 7) + f/2] ■ ( ' 

In the above, we have taken a uniform limit for the energy partition: p(z) = 1. 

Here, we have defined (e e ) = (E e )/N e = mean electron energy; {e p h) = {E p h) /N p h = mean phonon energy 
(= fesTg); and p e - P h = N e /N p h = electron-to- phonon number ratio. It is to be noted that in the limit a — 1 (i.e., 
for elastic electron-electron collisions as is usually expected for an electronic system unlike the case of the granular 
gas), and g ex — {i.e., no photo-excitation), we recover (e e ) = (e p /j), i.e., the electrons and the phonons are at the 
same temperature, as is physically expected under equilibrium conditions. In general, however, the mean electron 
energy in the steady state is not the same as the mean phonon energy, and the former depends on the excitation rate 
(the drive g ex )- 

The phenomenological integro-differential equation (Eq.38) is a general result. Note that in (Eq.38), the probability 
p(z) can be any 'collision distribution' function. This kinetic equation can be solved by numerical simulations with 
specific choice of parameters. The important aspect of this stochastic model is that it describes the kinetics in 
non-degenerate electronic sub-system, because in this regime microscopic calculation for non-equilibrium distribution 
functions are complex, and no simple model is available. Our simple stochastic model gives a partial solution to the 
problem. 



11 



V. CONCLUSION 



The present paper reviews the Two- Temperature Model based on Boltzmann-Bloch-Peierls kinetic equation. Its 
approximate validity is justified due to the fact that, in an excited metal, the intra-electron and intra-phonon relaxation 
time scales are much much less than the inter-electron-phonon relaxation time scales. In spite of the success of TTM 
in bulk case, it cannot be directly applied to spatially restricted systems, because of the non-resonant nature of 
the interaction between bulk phonons and the electrons. The main channel of energy relaxation is through electron 
surface-phonon interaction. Problems can be solved by considering surface effects, as explained in section III. We have 
seen (from the discussions in section IV) that TTM fails when one probe the sample with very short laser pulses (pulse 
width in the femto-second regime), and the non-thermal electron distribution has been experimentally observed. In 
this regard, we explained the rate-equation model of Sun etal and presented a stochastic kinetic model to describe 
the kinetics of non-thermal part of the electron distribution. So we see that TTM is an approximation of BBP kinetic 
equation. In fact, BBP approach also has its own limited validity, it is valid for a rarefied electron gas and in the 
first Born approximation. Its success here lies in the fact that conduction electrons in metals can be treated free. For 
a general non-equilibrium electron-phonon system, non-equilibrium statistical operator method should give us more 
profound results, that will go beyond all the cases considered here. 
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